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In this paper we examine to what extent traces of a clustering can be found in mean-field ground 
states of na-nuclei from 8 Be through 36 Ar as well as in some superdeformed states in 32 S, 36 Ar, 
and 40 Ca. For this purpose we calculate the overlap of the mean-field Slater determinant with one 
containg pure Gaussians and perfect spin and isospin symmetry, optimizing the overlap by varying 
the a-particle positions and radii. In some cases a coherent sum over different configurations is also 
employed. We find quite large overlaps for some of the lighter systems that diminish for nuclei above 
20 Ne but again strong clustering in 36 Ar. 



PACS numbers: 21.60.Gx, 21.60.Jz 
I. INTRODUCTION 

For light na nuclei, i. e., nuclei composed of 2n protons 
and 2n neutrons, a description in terms of a-clusters has 
enjoyed considerable success over many years. In their 
simplest form, these correspond to Mean-Field with a 
restricted wave function made up of a-clusters: quadru- 
ples of particles with both spin and isospin orientations 
and with Gaussian wave functions in space centered at 
given positions. These positions are then varied to obtain 
an optimal many-body wave function. Later variations 
of this method also allow variation of the radii, deformed 
clusters, and exploit the computational efficiency of the 
approach to calculate projected states (see pj and refer- 
ences therein). 

In the full Hartree-Fock or Mean-Field approach, on 
the other hand, the wave functions can be quite com- 
plicated in their spatial behavior and indeed in many 
modern calculations are freely variable over a coordinate 
grid, so that typically thousands of values are varied to 
optimize the many-body wave function. Since the Mean- 



Field state is a pure Slater determinant, it does not con- 
tain correlations in the conventional sense, which seems 
to exclude correlated objects like a-particles. A closer 
examination reveals, however, that correlations come in 
through the mean field. The identical wave functions 
assumed for each set of four particles in an a-cluster 
model are present in the full Mean-Field solution to a cer- 
tain extent, because the eigenfunctions in the same mean 
field potential will be filled with four particles, the main 
cause for deviations from the identical wave functions be- 
ing spin-orbit coupling (destroying spin degeneracy) and 
Coulomb effects (destroying isospin degeneracy). 
In this paper we address the following questions: 

1. To what extent do the Mean-Field states agree 
with an nev-cluster model, i. e. can they be rep- 
resented by cluster model wave functions? 

2. Does this agreement depend in a systematic way 
on spin orbit strength and the mass of the nuclei? 

3. Is the Mean-Field state better represented by a cor- 
related cluster state, in the sense of a coherent sum 
over different cluster configurations? 
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4. Are there indications for more exotic cluster states 
producible in the Mean-Field approach? 
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It is, of course, interesting to compare the results with 
cluster-model interpretations of these nuclei. In the fol- 
lowing, we will refer to the geometric configurations sug- 
gested by Zhang, Rae, and Merchant [2| based on a 
cranked version of the Block-Brink a-cluster model. In 
some cases also a comparison to deformed-basis antisym- 
metrized molecular dynamics (AMD) 0] calculations was 
performed. 

II. THE MEAN-FIELD CALCULATIONS 

A. Method of Calculation 

We represent the single-particle wave functions on a 
Cartesian grid with a grid spacing of 1 fm. The grid size 
is typically 24 3 for ground states and 36 x 24 2 for superde- 
formed states. This accuracy was seen to be sufficient to 
provide converged configurations. The numerical proce- 
dure is the damped-gradient iteration method [j], and 
all derivatives are calculated using the Fourier transform 
method. 

In the case of shape isomeric states, convergence is 
sometimes difficult to establish. Experience has shown 
that the observation of the change in total energy alone 
is not sufficient to judge convergence: this value usu- 
ally typically decreases rapidly to 8E/E k, 10~ 9 , which 
may indicate simply that the calculation is getting stuck 
and convergence is really stalled. A much better crite- 
rion turned out to be the mean-square deviations in the 
single-particle energies, summed over all states. These 
uncertainties thus measure how far the states still are 
from true eigenstates of the single-particle Hamiltonian. 

In practice it was often seen that the total energy 
seemed to converge perfectly well, while the fluctuations 
in the Hamiltonian stay at a relatively high level of 10~ 3 . 
In such cases usually the system will stay in this configu- 
ration for several thousand iterations and then transit to 
the ground state relatively rapidly. The underlying situ- 
ation can be visualized as a steep descent which arrives 
at some saddle point, where the direction of the descent 
must be modified strongly, keeping the system at this 
point for a very large number of iterations. Alternatively 
it may be caught in a shallow local minimum from which 
the accumulation of small changes in the wave functions 
enable the system to escape during further iterations. 

B. Interactions and symmetries 

We take four different Skyrme forces which all perform 
very well concerning nuclear bulk properties but differ in 
details: SkM* as a widely used traditional standard 
Sly6 as a recent fit which includes information on isotopic 
trends and neutron matter [(| , and SkI3 as well as SkM as 
recent fits which map the relativistic iso- vector structure 
of the spin-orbit force pj. SkI3 contains a fixed isovec- 
tor part analogous to the relativistic mean-field model, 



whereas SkM is adjusted allowing free variation of the 
isovector spin-orbit term. Thus all forces differ somewhat 
in their actual shell structure. Besides the effective mass, 
the bulk parameters (equilibrium energy and density, in- 
compressibility, symmetry energy) are comparable. In 
addition we consider two relativistic mean-field models. 
NL3 is a commonly used parameter set that includes 
non-linear interactions of the a mesons that generate the 
scalar interaction between the nucleons The model 
Xm is based on a chiral flavor-SU(3) effective Lagrangian, 
which is discussed in detail in 0. 

It will be seen that the results are generally in quite 
good agreement, with mainly the spin-orbit strength pro- 
ducing differences among the Skyrme forces; SkM* sys- 
tematically deviates somewhat from the others. For NL3 
and Xm the spin-orbit interactions are a direct conse- 
quence of the relativistic description of the nucleon wave- 
functions. The treatment of the center-of-mass motion 
is considerably different in SkM* from the other forces; 
since we do not include a center-of-mass correction, there 
are considerable differences in the binding energies of 
light nuclei. 

The wave functions were not subject to any symmetry 
requirements, having neither good spin nor good parity 
quantum numbers; in fact, the initialization purposely 
avoided symmetries to make sure that these were not 
accidentally maintained. 

The calculations included pairing in the BCS approxi- 
mation with a delta force, taking into account three times 
the number of occupied states without pairing. The 
parameters were the ones usually associated with each 
Skyrme force. 



C. Observables 

The principal observables that will be quoted are the 
total binding or excitation energy, the spin-orbit energy, 
the deformation, and the intrinsic quadrupole moment. 
The quadrupole moment was calculated in the standard 
way as 

Q = J d 3 rp{r){3z 2 -r 2 ) (1) 

for the case that the z-axis is the axis of axial symmetry. 
A deformation parameter is then obtained as 

A being the mass number and R the mean-quare radius 
of the system. This prescription is, of course, not unique, 
since the dependence on the mean-square radius is some- 
what arbitrary, but follows [l(j . Oblate deformations are 
recognizable by the negative sign in Q20 and /?20- 
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III. METHODS OF ANALYSIS 

To investigate possible molecular structure, three dif- 
ferent methods of analyzing the mean-field Slater deter- 
minant |\&) were employed. In this section the procedures 
are described, while detailed results for the various nuclei 
will be given in Section IVI and IVIl 



A. na-cluster model 

We investigated the clustering structure by construct- 
ing a model many-body wave function with clustering as- 
pects built in and maximizing its overlap with the Slater 
determinant obtained from the Mean-Field calculation. 
As the single-particle wave functions have no good quan- 
tum numbers aside from the energy, only the overlap of 
the complete Slater determinants can be meaningfully in- 
vestigated. In case of the relavistic models we used the 
upper components of the wavefunctions to calculate the 
overlap. 

In the simplest case, which we shall refer to as the 
na-model, the model wave function is constructed com- 
pletely out of a-particle configurations. For a nucleus with 
A nucleons, we use N a = Aj 4 clusters located at ffc, 
k = 1, . . . N a and with radius parameters Ofc. Depending 
on the total number of a-particles, we usually constrained 
their radii to be equal, but sometimes, when total com- 
puting time allowed, also allowed them to vary to see 
whether the overlap is thereby improved. The resulting 
improvement was of the order of a few percent typically, 
so that in most cases we will quote results only with one 
common value of a adjusted for maximum total overlap. 

The spatial part of the wave functions then is given by 



4> k (r) = A k exp 



(r- f fe ) 2 



2*2 



1 . 



■ N a . 



(3) 



with Ak a normalization factor that is computed numer- 
ically. The numerical normalization factor calculated on 
the same spatial grid with 1 fm grid spacing as the Mean- 
Field wave functions was found to agree with the analytic 
one to five significant digits, so that the spatial grid ap- 
pears sufficiently accurate for the overlap calculations. 
Each of these wave functions occurs four times in the 
model Slater determinant: for proton and neutron and 
with spin up or down. 

Note that the na-cluster single-particle wave functions 
are not orthogonal. As a consequence, the norm of the 
model Slater determinant must be calculated as the de- 
terminant of the single-particle overlaps 



A/- = det«0*|0 fc }) 



(4) 



This determinant can become very small if the wave 
functions have large overlaps, which occur frequently in 
the optimal configurations, indicating that the cluster 
description does not correspond to clearly separated a- 
particles, but generates the Mean-Field states largely by 



antisymmetrization. For example, two Gaussians with 
the centers placed very close to each other generate an 
odd-parity state through antisymmetrization. 

The overlap O = |(^|3>) | 2 given in the tables below 
is always the absolute square of the matrix element be- 
tween the Mean-Field Slater determinant and the nor- 
malized na-cluster wave function |$), so it describes the 
probability for the Mean-Field state to agree with the 
cluster model state. Since the Mean-Field wave func- 
tions have good isospin but are mixed in spin projection, 
the total overlap could be simplified by decomposition 
into a proton and a neutron factor. Although the effects 
of Coulomb are practically negligible, the calculation did 
treat the two isospins separately. 



B. Calculation of na-cluster configurations 

The most effective method to search for the optimal 
positions and sizes of the na-cluster positions was a 
search with an optimization algorithm starting from ran- 
dom positions (within a box roughly corresponding to the 
nuclear volume) and a standard starting value for the ra- 
dius parameters of 1.8 fm. We found that optimizing the 
logarithm of the overlap of the two Slater determinants 
with the algorithm of Stewart ^l| (available as routine 
"dmnfb" in NETLIB) led to very rapid convergence even 
when the individual particles were placed very badly ini- 
tially with a starting overlap of < 10 -15 with the Mean- 
Field Slater determinant. 

The final configurations could in all cases be decom- 
posed into 

1. isolated a-particles in a very well determined posi- 
tion characterizing the cluster geometry, and 

2. groups of a particles spaced very closely at dis- 
tances typically < 0.1 fm. These represent cores 
made up of larger nuclei, usually 12 C or 16 0. Since 
the wave functions in these clusters serve mostly 
to generate higher shell-model states through an- 
tisymmetrization, their geometrical layout appears 
random. It was found that the overlap with the 
mean-field wave functions changes only by about 
2% when the distance between the particles in- 
creases from 0.01 fm to 0.8 fm. It is largest for 
the smallest distance, but this can cause numeri- 
cal problems because of the nearly degenerate wave 
functions. 

It should be noted that the positions found usually 
showed some alignment with the Cartesian axes even for 
spherically symmetric nuclei, since the representation on 
a Cartesian grid violates that symmetry to a small but 
noticeable degree. 

In some cases a fixed geometry was also employed to 
check how different cluster geometries compete in de- 
scribing the HF ground state. In these cases a small 
number of geometric quantities was varied to achieve op- 
timal overlap. 
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C. Model of core nucleus with additional 

a-particles 

In some cases, a model wave function was used where 
a core of 12 C or 16 replaced the corresponding number 
of a-particles. Note that the wave functions used for 
these were simply static Mean-Field solutions for those 
core nuclei, without any adjustable parameters (except 
for the position of their centers). In the cases where 
three or four a-particles tended to contract into the same 
location, this may produce more stable numerical results 
because of the norm going to zero for the a-particle case. 
The most interesting application of this technique is to 
the interpretation of the strongly deformed state of 32 S 
as an 16 0+ 16 configuration. 

D. Overlap with a collective cluster space 

In some cases it was tried to allow some parameter 
determining the na-cluster configuration to vary in the 
model wave function and to calculate the overlap with 
the collective space generated in this way. As an exam- 
ple, in the case of the strongly deformed 32 S with Slater 
determinant state interpreted as an 16 0+ 16 con- 
figuration, the Slater determinants consisting of the 
antisymmetrized combination of two 16 O nuclei at dis- 
tances di, i = 1 . . . TV were used as a nonorthogonal basis. 
The matrix elements of the overlap kernel for this basis 
are A i} = (<J>j|$j). The probability for the state \&) to 
be within this space is given by 

N 

°= E (^-%<*il*)(*|3i>. (5) 

The result depends on TV, but converges rapidly as soon 
as the collective space is sufficiently spanned. The value 
of TV of course needs to be chosen carefully. We found 
that numerically the results remain stable for TV < 20, 
but usually n ~ 5 is already sufficient for convergence. 
The criterion of Eq. J5J to extend the Slater stgate 
can be reproduced by a coherent superposition of 
a-model states 

|^(modcl) ) = J2 |$ 2 ) Cj (6) 

i 

with optimized coefficients Cj. This ansatz is the discrete 
version of the celebrated Generator Coordinate Method 
(GCM) ^3] which had, in fact, one of its early success- 
ful applications in the realm of the a-cluster model |l3j . 
Note that the criterion © compares somewhat different 
concepts. The Mean-Field state \$f) is the optimum one 
can achieve in the space of pure Slater states. The GCM 
ansatz on the other hand, is restricted to Slater states 
within the a-model but goes beyond Slater space in that 
it can include substantial correlations in its coherent su- 
perposition. In spite, or just because, of this conceputal 
difference, it is most interesting to investigate the mutual 
overlap of these two lines of approach in Eq. (J5J • 



IV. GENERAL FEATURES OF THE RESULTS 

A general feature of all nuclei considered was that the 
ground states are axially symmetric as well as reflection 
symmetric for all the forces considered, although nei- 
ther the initialization nor the computation enforced these 
symmetries.. 

We found in all cases that the converged station- 
ary states showed no pairing, so that calculations with- 
out pairing might appear sufficient. There is, however, 
the danger that in iterations using the restricted single- 
particle space of the occupied wave functions the two 
fragments may get stuck in a local configuration. It is 
known that without pairing the deformation energy curve 
is an envelope of many relatively steep curves correspond- 
ing to specific configurations and not all of the narrow 
minima produced for such a configuration remain min- 
ima in the envelope. 

In practice it was found that only part of the highly 
deformed states were still produced in the calculation if 
pairing was included. Although again in the final result 
pairing disappeared, we therefore discuss only these ulti- 
mately stable configurations. 

The freedom allowed by the three-dimensional code 
was found not to be necessary in these pure Mean-Field 
cases, as the stationary states in all cases turned out to 
axially symmetric. A comparison with na-cluster model 
wave functions of course required three-dimensional ge- 
ometry. A code with axial symmetry and the facility for 
quadrupole constraints was used in some cases to support 
the results and check deformation dependence. 



V. GROUND STATES OF na NUCLEI 

In the following we summarize the properties of the na- 
nuclei below 40 Ca. Density contour plots are shown only 
for cases with interesting structure and only for SkI3; in 
these the contour lines are spaced by 0.01 fm -3 . The con- 
tour lines sometimes appear non-smooth; this is because 
the graphics program does not interpolate. The physical 
wave functions are quite smooth owing to Fourier expan- 
sion. In all tables, we use Ei s for the spin-orbit energy, 
and O for the overlap with na- or collective cluster con- 
figurations. 

For some specific nuclei we also mention the overlap be- 
tween the Mean-Field wave function and the AMD wave 
function, and also that between AMD and the na-cluster 
wave function. Here the AMD wave function is obtained 
using the Gogny D1S force 0. The use of a different 
type of force will, of course, reduce the agreement with 
the Mean Field calculations, but a comparison is never- 
theless instructive. 
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TABLE I: Binding energies in MeV calculated for the nu- 
clei considered in this paper with the four Skyrme force 
parametrizations used. 
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27.8 


49.7 
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SkI4 


27.7 
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128.6 
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155.9 
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Exp. 


198.3 


236.5 


271.8 


306.7 




SkI3 


194.7 


233.0 


267.1 


304.8 




SkM 


196.11 


234.9 


269.4 


305.5 




SLy6 


194.4 


233.1 


268.5 


304.0 




SkM* 


197.6 


237.9 


275.1 


305.5 




NL3 


194.1 


231.8 


265.6 


302.2 




Xm 


197.0 


235.0 


269.5 


308.8 





A. Binding Energies 

In Table [i] the binding energies obtained in all of the 
ground-state calculations are listed for reference together 
with the experimental values. All calculations contain 
a large contribution from the center-of-mass correction, 
which is calculated in different ways: the Skyrme forces 
SkI3, SkM, Sly6 and the % m -model contain a microscopi- 
cally calculated correction, NL3 uses the harmonic oscil- 
lator approximation, and SkM* a correction to the nu- 
cleonic mass m — > m — m/A. 



B. Cluster configuration overview 

The optimal cluster configurations obtained are de- 
scribed below in detail. To facilitate visualization, we 
show an overview of the geometry in Fig.^ Those cases 
where there is a trivial configuration as well as those with- 
out clear structure preference are omitted. 

C. The Nucleus 4 He 

The a-particle itself is of course an extreme case for a 
mean-field description. We summarize its properties in 
Table im mostly because it allows to judge how justified it 
is to use Gaussian wave functions for the overlap calcula- 
tions with an na-cluster structure in heavier systems and 
how the optimum radius depends on the force. For calcu- 
lating the overlap of the a-particle itself with the Gaus- 
sian model wave functions, only the radius parameter a 
was varied with the optimum value given in the table. In 
all cases the overlap is remarkably close to one, showing 
that the use of the Gaussian is a good approximation. 
Note that in the case of the relativistic descriptions the 
width of the Gaussian is somewhat smaller than for the 




FIG. 1: Visualization of the cluster geometry present in the 
Mean-Field ground or excited superdeformed states (denoted 
by asterisks) for the nontrivial cases. The figures only illus- 
trate the topology; no attempt was made to reproduce correct 
scale lengths. Double links between particles indicate very 
small distances. The triangles and pentagon are perpendicu- 
lar to the linear links. 

TABLE II: Physical properties of 4 He for the different forces 
considered together with data describing the optimum repro- 
duction by Gaussian wave functions. Note that for the rela- 
tivistic cases there is no separate LS-force, so no Ei 3 values 
are given. 



Force 


Ei s 


Rrvns 


Gaussian Model 




[MeV] 


[fm] 


O [%] 


a [fm] 


SkI3 





2.06 


99.5% 


1.67 


SkM 





2.07 


98.8% 


1.68 


Sly6 





2.12 


98.8% 


1.72 


SkM* 





1.97 


99.4% 


1.61 


NL3 




1.87 


99.4% 


1.54 


Xm 




1.89 


89.8% 


1.54 



Skyrme forces, which leads to smaller widths throughout 
the a-cluster analysis presented in this section. 

D. The Nucleus 8 Be 

This nucleus is a clear-cut case for a molecular inter- 
pretation. For all Skyrme forces considered it turned out 
to have a necked-in shape agreeing well with a double-a 
structure interpretation. Density isolines illustrate this 
in Fig. [3 while other properties are summarized in Ta- 
ble IIIII The physical properties are quite similar, with 
SkM* being the odd force out, as expected. 

The na-cluster analysis naturally found two a- 
particles positioned symmetrically along the symmetry 
axis with a distance d and a radius of a. The results 
in the table indicate that this interpretation works ex- 
ceedingly well, as even the radius parameters a agree re- 
markably well with those for the free a-particle. There is, 
however, a notable nonorthogonality of the two clusters. 
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FIG. 2: Density contour lines for Be. 



TABLE III: Physical properties of Be for the four Skyrme 
forces considered. For the cluster model interpretation the 
overlap is given together with the optimum distance and ra- 
dius of the clusters. The last two columns indicate the overlap 
with an q-GCM model state and the optimal a obtained. 



Force 


Eis 


02 


Q20 


Cluster Analysis 


q-GCM 




[MeV] 




[fm 2 ] 


0[%] 


d[fm] 


[fm] 


0[%] 


a [fm] 


SkI3 


2.5 


0.677 


46.8 


82 


2.70 


1.68 


98 


1.65 


SkI4 


2.7 


0.669 


46.2 


81 


2.64 


1.69 


97 


1.64 


Sly6 


2.4 


0.666 


48.0 


81 


2.68 


1.73 


97 


1.68 


SkM* 


4.6 


0.593 


38.0 


71 


2.20 


1.68 


82 


1.62 


NL3 




0.679 


40.7 


85 


2.52 


1.56 










0.671 


41.8 


68 


2.60 


1.60 







This nucleus also gave a good opportunity to test the 
overlap with a-GCM cluster wave functions as described 
in section ITlIDI We first used the distance d as a collec- 
tive coordinate, distributing 10 points in the interval be- 
tween d — 1 fm and d = 5 fm and varying a to maximize 
the overlap. The result is given in the "a-GCM" column 
in the tables. The overlap increases noticeably. This 
calculation was performed only for the Skyrme- forces. 

It is remarkable that the optimum a values are smaller 
for the "a-GCM" case. It may be that for the single- 
distance analysis the a-particles have to widen a bit to 
simulate the effect of orbital motion; one should not, how- 
ever, read too much into this as the difference in overlap 
between the different a values is of the order of 3%. 

The AMD wave function also has a quite large over- 
lap of 96% with the collective two-a-cluster model space. 
The overlap between the AMD and the Mean-Field wave 
function is a little bit smaller at 85%. 



E. The Nucleus 12 C 

The surprising feature in this nucleus is the discrep- 
ancy in ground-state properties. While three of the 
Skyrme forces used have a spherical ground state, SkI3 
produces a strongly deformed oblate shape. In order to 
understand this, a constrained axially symmetric calcula- 
tion was also performed, which showed that the structure 
is really quite similar in all cases: the potential is very 
flat as a function of deformation, and it just happens 
that the oblate point is about half an MeV lower than 
the spherical one for SkI3, while this is reversed for the 



other three forces. 

TABLE IV: Physical properties of 12 C for the four Skyrme 
forces considered. The other columns give the maximum over- 
lap with an a-particle configuration for a fixed triangle with 
side length 0.01 fm. 



Force 




02 


Q20 


Cluster Analysis 




[MeV] 




[fm 2 ] 


0[%] 


a 


SkI3 


15.8 


-0.256 


-24.7 


28 


1.72 


SkM 


21.6 


0.000 


0.00 


1 


1.66 


Sly6 
SkM* 


19.8 


0.000 


0.00 


1 


1.68 


23.6 


0.000 


0.00 


1 


1.64 


NL3 




0.000 


0.00 


1 


1.48 


Xm 




0.129 


11.4 


2 


1.64 



The search for the optimal na-cluster configuration in 
all cases yielded the a-particles very close to the center of 
mass of the nucleus in a roughly triangular configuration. 
The overlap is much larger for the oblately deformed nu- 
cleus, since there is a preferred spatial alignment for the 
triangle in this case. A calculation with fixed triangular 
geometry yielded very similar results. 

An interpretation with an a-GCM cluster configura- 
tion was attempted in two ways. First, a number of dif- 
ferent sizes of the triangles were used. The configuration 
which gave converged results was 10 points distributed 
between 0.01 fm and 3 fm. This improved the overlap 
very little, for example to 30% for SkI3. 

Alternatively rotated configurations were included. 
Adding the Slater determinants for the triangle rotated 
in its plane by 2-10 equally spaced angles changed the 
overlap by less than 1%, adding three-dimensional rota- 
tions also produced similarly small effects. 

The overall disappointingly small overlap with a- 
clusters in this nucleus appears to be due to the large 
spin-orbit energy, with only the deformed state for SkI3 
a better candidate for clustering. Because the oblate and 
spherical states are close in energy and not separated by 
a barrier, including vibrational correlations in the mean- 
field model should lead to an intermediate value for the 
overlap that is quite comparable for all forces. 

F. The Nucleus 16 O 

For this nucleus all forces produce a spherical ground 
state, so that in this case the mean-square radius is also 
an interesting quantity to compare. It is shown in Ta- 
ble and it is remarkable that the radii agree very well 
for all four forces. 

The cluster analysis with free placement produced sim- 
ilar results to the spherical 12 C nuclei: the a-particles 
tended to the same position at the center off mass, indi- 
cating that one is not dealing with well-separated clusters 
but again with wave functions produced mainly by anti- 
symmetrization. The optimization located the particles 
at positions (d, 0, 0), (0, d, 0), (0,0, d), and (— d, 0, — d) 
with d — 0.01 (as mentioned, results are insensitive to 
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TABLE V: Properties of 1 for the four Skyrme forces to- 
gether with results for free cluster placement. 



Force 


Eis 


Rrms 


Cluster Analysis 




[MeV] 


[fm] 


0[%] 


cr [fm] 


SkI3 


1.0 


2.65 


96 


1.76 


SkI4 


1.0 


2.65 


96 


1.76 


Sly6 


0.9 


2.69 


96 


1.79 


SkM* 


1.1 


2.68 


96 


1.78 


NL3 




2.56 


95 


1.72 






2.58 


79 


1.71 



the exact value of d). This is not a natural arrangement, 
but seems to indicate that the principal requirement is 
that the four particles not be coplanar. The small varia- 
tion in the tr-values reflects that in the mean-square radii. 

A cluster analysis was attempted using a tetrahe- 
dron of a-particles, with the particles' radii and mu- 
tual distance as parameters. The results was that max- 
imizing the overlap drives this distance towards zero at 
(j w 1.76fm. In this case the limiting overlap is quite close 
to 1, showing that the tetrahedral symmetry is a good 
basis for expanding the true wave functions. Like in the 
case of 12 C, however, the cluster interpretation is thus 
invalid in this case. It should be noted that the overlap, 
which reaches 98% for a distance of 0.01 fm, only goes 
down to 97.5% at 0.2 fm. This opens the alternative of 
either using 16 O wave functions or the four-a structure 
to look for the presence of an 16 O core in heavier nuclei. 

That the precise arrangement of the clusters is not im- 
portant was also demonstrated by using fixed geometry. 
Using an exact tetrahedron configuration produced iden- 
tical results; only when the particles were coplanar the 
overlap vanished immediately. 

Our results are in complete agreement with the cluster 
model of @|. 

In view of the good description of this nucleus by static 
cluster wave functions we did not perform multiconfigu- 
ration studies in this case. 




1 3 5 

x [fm] 



FIG. 3: Density contour lines for Ne. 



G. The Nucleus 20 Ne 

In this case the shape of the nucleus, shown in Fig. [3] 
shows a strong prolate deformation with a sugestion of a 
central, more or less spherical part, with two additional 
mass distributions added on both sides. The physical 



properties are given in Table IVII This would suggest a 
12 C core with two a-particles added, and this is indeed 
what was produced by the unconstrained search. This 
is quite close to the interpretation in 0, although there 
the central 12 C has an annular structure. 



TABLE VI: Physical properties of Ne for the four Skyrme 
forces considered. The parameters for the cluster analysis are 
described in the text. 



Force 


Eis 




Q20 




Cluster Analysis 




[MeV] 




[fm 2 ] 


0[%] 


d [fm] 


<ri [fm] 


02 [fm] 


SkI3 


8.5 


0.423 


91.0 


53 


1.91 


1.78 


1.71 


SkM 


9.2 


0.412 


88.0 


49 


1.86 


1.78 


1.71 


Sly6 
SkM* 


8.5 


0.409 


89.8 


47 


1.84 


1.80 


1.74 


11.1 


0.371 


79.2 


36 


1.59 


1.77 


1.73 


NL3 




0.425 


84.4 


59 


1.83 


1.74 


1.65 






0.439 


90.6 


50 


1.89 


1.74 


1.69 



It appeared expedient to try different radii for the a- 
particles in this case. Letting the program optimize all 
the radii simultaneously instead of using a common a 
led to an improvement of < 2% in overlap and naturally 
showed the three particles in the core with the same ra- 
dius, while the two displaced ones also were identical. 
The final configuration, for which details are listed in the 
table, thus has three a-particles with radius <j\ in a tri- 
angle at the central position and perpendicular to the 
longest axis of the nucleus, while two with radius 02 are 
positioned at id along that axis. 

An attempt with a multiconfiguration analysis yielded 
disappointing results: rotating the central triangle gave 
less than 1% improvement, while allowing the outer par- 
ticles to move along the axis yielded up to 2%. Adding 
different sizes of the central triangle gave similar values. 

The same 5a-configuration is also obtained in an AMD 
calculation without parity projection before variation. 
The overlap with the pure na configuration amounts to 
57%. A larger overlap of 61% is obtained by including 
parity projection before variation. In this case, the op- 
timum configuration corresponds to an a+ 16 cluster 
nature. The overlap between the AMD and Mean-Field 
wave functions is 78% and 37% without and with parity 
projection, respectively. 



H. The Nucleus 24 Mg 

This nucleus is characterized by a strong quadrupole 
deformation for all forces considered, as seen in Tabic lyTTI 
and in Fig. ^ 

The unconstrained positioning of the a-particles led to 
a grouping in doublets in the following way: with the x- 
axis being the symmetry axis, two doublets are arranged 
at x = ±p with the particles displaced at slightly positive 
and negative y-values (around ±0.01), while the third 
doublet is obliquely shifted from the axis to ±(0, q, r) 
with p and q surprisingly at larger values. While this 
appears to show some indications of true clustering, the 
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overlap is disappointingly small. The oblique positioning 
of the central doublet again indicates the need for having 
the spatial degrees of freedom sufficiently represented. 



TABLE VII: Physical properties of 24 Mg. The parameters for 
the cluster analysis are described in the text. 



Force 




fa 


Q20 




Cluster Analysis 






[MeV] 




[fm 2 ] 


0[%] 


P [H 


q [fm] 


r [fm] 


a [fm] 


SkI3 


22.6 


0.423 


117.2 


2.3 


1.40 


0.17 


0.37 


1.77 


SkI4 


23.5 


0.416 


114.1 


2.2 


1.38 


0.10 


0.41 


1.76 


Sly6 


21.6 


0.413 


116.5 


2.2 


1.39 


0.04 


0.40 


1.79 


SkM* 


25.4 


0.389 


107.0 


1.8 


1.31 


0.15 


0.41 


1.77 


NL3 




0.416 


103.9 


2.0 


1.37 


0.17 


0.45 


1.71 






0.431 


115.2 


2.0 


1.38 


0.25 


0.31 


1.74 



the situation turns around: the overlap with the pure na- 
clustcr wave function is reduced to a few percent, while 
that with the Mean-Field one goes up to 68%. 



I. The Nucleus 28 Si 

Except for SkM*, which shows a relatively pure 
quadrupole deformation, all forces have this nucleus 
oblate deformed but with a strong hexadecupole con- 
tribution visible in the high-density contour lines (see 
Fig: |3J. The oblate deformation is in agreement with 
other calculations [Tol IT1 ITl \l% . 

As suggested by this geometrical shape, the optimal 
placement resulted in 5 a-particles in a pentagon per- 
pendicular to the symmetry (z-) axis and two further 
particles at symmetric positions on either side. The dis- 
tances between the particles again tended towards zero, 
and because of that the pentagon was not a regular poly- 
gon. Doing the calculation with a regular pentagon gave 
identical results for the overlap. The Mean-Field wave 
functions thus seem to indicate a preference for the pen- 
tagon symmetry, especially in view of the fact that the 
overlaps are quite large in this case. 

TABLE VIII: Physical properties of 28 Si and optimum cluster 
parameters. 



Force 


Els 


P2 


Q20 


Cluster Analysis 




[MeV] 




[fm 2 ] 


0[%] 


a [fm] 


SkI3 


30.2 


-0.318 


-108.5 


7.8 


1.86 


SkU 


33.6 


-0.294 


-98.6 


4.5 


1.85 


Sly6 


31.7 


-0.289 


-98.8 


4.0 


1.87 


SkM* 


40.6 


-0.224 


-73.7 


0.8 


1.83 


NL3 




-0.303 


-93.5 


6.2 


1.80 






-0.351 


-117.9 


14.9 


1.85 



TABLE IX: Physical properties of 3 S for the four Skyrme 
forces considered. 



Force 


Eis 


fh 


Q20 


Cluster Analysis 




[MeV] 




[fm 2 ] 


0[%] 


a 


SkI3 


29.6 


0.222 


90.6 


1.6 


1.87 


SkU 


37.0 


0.147 


58.7 


0.3 


1.86 


Sly6 


40.0 


0.000 


0.0 


0.04 


1.87 


SkM* 


44.7 


0.000 


0.0 


0.03 


1.86 


NL3 




0.228 


86.24 


1.8 


1.83 






0.253 


101.1 


2.2 


1.85 



J. The Nucleus J2 S 

This nucleus again shows a strong dependence on the 
force. For SkI3 and SklA it is prolate but with the strange 
feature of an oblate central part at high densities (see 
Fig. |SJ), whereas it is spherical for the other two forces. 
Other self-consistent calculations fp | .[Tol 111 HI lv \ ) 





FIG. 4: Density contour lines for 24 Mg (left) and 28 Si(right). 



Alternative interpretations with fixed geometrical 
shape were also explored: 



1. 



2. 



3. 



Two equilateral triangles displaced to x — id and 
either aligned the same way perpendicular to the 
symmetry axis or rotated by 30° relative to each 
other. In both cases the maximum overlap that 
could be achieved by varying the triangle size and 
d was < 0.016. Taking the overlap with two 12 C 
nuclei displaced along the symmetry axis yielded 
very small overlap (< 10~ 9 ). 

Placing a square with vertices at x = 0, y, z = 
±g and two additional particles at x — ip yielded 
microscopic overlaps of about 0.0002. Clearly this 
arrangement does not catch the symmetry of the 
single-particle space near the center. 

Attempts of putting an 16 nucleus into the center, 
accompanied by two a-particles on both sides, lead 
to much smaller overlap. 



The interpretation of Ref. [2( in this case does not agree 
well; they favor the second of the three choices given. The 
large value of the spin-orbit energy in our case seems to 
make a cluster expansion less applicable. 

The strong effect of the spin-orbit force to dissolve the 
cluster structure is also confirmed by the AMD analysis. 
When the spin-orbit force is switched off, the AMD wave 
function has a 2a+ 16 cluster structure and large over- 
lap (56%) with the na-cluster wave function, but a small 
overlap of a few percent with the Mean-Field wave func- 
tion. When the spin-orbit is included, on the other hand, 



9 




x [fm] x [fm] 



FIG. 5: Density contour lines for 3 S and Ar. 



generally prefer a prolate deformation. The uncon- 
strained search for an a-configuration in this case pro- 
duced very small overlaps; the geometric positions in all 
cases appeared almost random, filling a cube near the 
center of mass. The distances, however, in this case were 
up to the order of 0.5 fm, showing that the overlap is 
quite insensitive to the particle placement. The very 
small overlap may be due to the lsi/2 filling and the 
large contribution of the spin-orbit force. 

The AMD wave function in this case is closer to the 
Mean-Field one with an overlap of 40%, but has only a 
few percent overlap with the na-cluster wave function. 



K. The Nucleus 36 Ar 

This nucleus appears as just the opposite of 32 S: it is 
oblate with a deformation in agreement with other cal- 
culations 0, 0, 0, 0| , but shows an inner core of pro- 
late deformation at high density. The geometry of the 
na-cluster configurations found appeared as random as 
in the case of 32 S, but the distances to the center-of-mass 
were again of the the order of < 0.1 fm, which together 
with the much larger overlaps shows that is again is a case 
of a dominating antisymmetrization effect. This comes 
close to the almost spherical geometry preferred by 

It is, of course, somewhat unexpected that in this nu- 
cleus clustering appears stronger than for all the others 
heavier than 20 Ne. This should be due to the much re- 
duced role of the spin-orbit coupling in this nucleus as 
witnessed by the spin-orbit energy being only half of that 
for 32 S. 



TABLE X: Physical properties of Ar and cluster analysis. 



Force 


Eis 


#> 


Q20 


Cluster Analysis 




[MeV] 




[fm 2 ] 


£>[%] 


a 


SkI3 


13.5 


-0.183 


-90.3 


28 


1.92 


SkI4 


15.2 


-0.171 


-83.8 


21 


1.91 


Sly6 


15.1 


-0.162 


-81.0 


17 


1.93 


SkM* 


19.1 


-0.133 


-65.9 


7 


1.92 


NL3 




-0.186 


-87.3 


33 


1.90 


Xm 




-0.194 


-94.1 


29 


1.89 



VI. STRONGLY DEFORMED STATES 
A. Technique for searching such states 

The search for strongly deformed isomeric states in the 
Mean-Field approximation has usually been done with 
constraints. We tried a new method that has the po- 
tential for producing more exotic configurations: initial- 
ization of the Mean-Field calculation with a fragment 
configuration. Thus, e. g., one may place two 16 nuclei 
at a certain distance d from each other into a numerical 
grid, orthogonalize, and then iterate to converge to a sta- 
tionary state. For small distances it may converge into 
the ground state — although we find this rarely happens, 
because the antisymmetrization produces quite a differ- 
ent set of single-particle orbitals than that of the ground 
state — , for a large d the nuclei may be driven apart, 
and finally, for some range convergence may happen into 
a highly deformed configuration. 

The attractive feature of this method is that initial 
configurations are highly flexible and might lead into 
configurations that cannot easily be reached with a con- 
straint. For example, a chain or a triangle of three a- 
particles may be used to look for a strongly deformed 
exotic state in 12 C, or one might try various asymmet- 
ric mass configurations to to investigate whether they all 
converge to the same deformed state. 

With pairing omitted, we found a large number of ex- 
otic states, which, however, corresponded to local minima 
under the constraint of a fixed single-particle configura- 
tion. Such minima may be important for the dynamics, 
as they correspond to large collective mass when the col- 
lective motion in the deformation is slow enough to cross 
over to another single-particle configuration, but clearly 
are not isomeric states. We therefore used pairing in the 
iterations, which made most of the minima disappear; 
the final states, however, as for the ground states, did 
not show significant pairing. 

B. General features of superdeformed states 

A large number of different combinations of na-nuclei 
leading to compound systems up to Ca were considered 
and the static calculations with different initial fragment 
positions were performed. Before going into the detailed 
results given in Table IXT1 and the following sections, let 
us summarize the main points of the findings: 

• Superdeformed shapes were found in the majority 
of cases, but none of these appears to be of clearly 
molecular type. 

• These exotic states were always axially and 
reflection-symmetric, even if the calculation was 
started with two nuclei of different mass. 

• For the superdeformed states only Skyrme forces 
SkI3 and Sly6 were employed. Generally the two 
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TABLE XI: Physical properties of highly deformed states for 
the Skyrme force SkI3 and the chiral model Xm- 





E* 


Eis 




Q 2 


Cluster Analysis 


System 


[MeV] 


[MeV] 




[fm 2 ] 


0[%] 


a [fm] 


SKI3 


28 Si 


13.2 


25.6 


0.773 


325.9 


1.3 


1.80 




8.9 


12.3 


0.737 


377.4 


34.5 


1.79 


36 Ar 


8.7 


29.1 


0.529 


283.4 


3 


1.88 


40 Ca 


26.3 


28.1 


0.983 


859.8 


0.8 


1.80 




28 Si 


11.1 




0.775 


314.4 


1.2 


1.77 


32 g 


5.7 




0.743 


368.6 


36.2 


1.75 


36 Ar 


9.2 




0.533 


278.9 


1.2 


1.84 


40 Ca 


23.7 




0.982 


829.41 


0.8 


1.77 



Skyrme forces yielded similar results, with the exci- 
tation energy somewhat more force-dependent than 
the deformation. We only present results for SkI3 
and mention discrepancies to Sly6 where they ap- 
pear interesting. 

• For these cases a comparison to the results of 
Ref. is not given: although the densities are 
roughly similar at least in the case of 32 S, they did 
not use a spin-orbit force, which makes the com- 
parison of the cluster interpretations problematic. 

The results are summarized in Table IXll 



C. The Compound System 24 Mg 

In the combination of two 12 C nuclei, no indication of 
a superdeformed state was found. The iterations always 
lead either to the ground state of 24 Mg or to further 
separation. 



D. The Compound System 28 Si 

This is formed in the combination 12 C+ 16 0. The it- 
erations lead - aside from the ground state of 28 Si- to 
one superdeformed state. Its properties, like of all the 
states found, are summarized in Tabic IXll The density 
contour lines show a bit of a double-humped structure 
at high density (see Fig. EJ). The unconstrained cluster 
analysis showed that along the symmetry axis the parti- 
cles are grouped in the following way: one is located at 
a central x-position but with a half-fermi displacement 
in the y and — z directions (this displacement is proba- 
bly accidental and not significant for the overlap, since it 
violates overall reflection symmetry). At ±1.23 fm there 
are pairs of particles on both sides of the axis at a rela- 
tive distance of about 0.1 fm, while finally two particles 
are located at ±3.1 fm on the axis. The overlap with the 
na-cluster configuration is, however, quite small. 




FIG. 6: Density contour lines for Si and Ssuperdeformed 
states. 



E. The Compound System 32 S 

This was studied in two different configurations: 
12 C± 20 Ne and 16 0+ 16 0, so that an interesting aspect is 
whether the initial asymmetry is lost. It was found that 
both configurations converge rapidly into the same state, 
which is shown in Fig. A state of similar deformation 
and excitation ene rgy is well known in the literature with 
similar properties 0, 120L l2ll I22I I23I ] . It shows a slight 
amount of necking-in, making this the best candidate for 
a "molecular" configuration. 

The relatively small spin-orbit energy, compared to 
the other superdeformed states, makes a cluster analy- 
sis much more attractive. The unconstrained cv-analysis 
shows all particles very close (the usual 0.01 fm) to the 
symmetry axis with two triplets at ±1.9 fm and addi- 
tional single particles at ±2.82 fm. The center-of-mass 
of the four particles on each side is at ±2.1 fm. This is a 
clear indication of a molecular structure with somewhat 
distorted 16 O nuclei at a distance of about 4 fm. 

For this reason an analysis of this system based on 
a double 16 configuration was also attempted. Two 
ground-state configurations were inserted into the grid at 
a specified distance, with the two nuclei being displaced 
symmetrically from the center-of-mass of the compound 
system along the axis of symmetry. The resulting over- 
lap is shown in Fig. \7\ together with the norm of the 
combined Slater determinant. The resulting maximum 
overlap indicates that this molecular configuration is not 
a perfect match; it is only near 6% and, moreover, falls off 
relatively slowly as the distance is varied. On the other 
hand, the norm of the model wave function (after or- 
thogonalization) is already close to unity, indicating the 
possibility of molecular structure at this distance, and 
using undistorted 16 O nuclei for the overlap analysis is a 
strong constraint. 

The small overlap, however, is not the full story. Do- 
ing an analysis along the lines of section IIII Dl with 20 
different c. m. distances used between 3 and 8 fm, we 
find a total overlap with this collective space of 35% - 
which agrees surprisingly well with the unconstrained a- 
particle analysis. The much more comprehensive anal- 
ysis of this system done by Kimura and Horiuchi [24| 
gives a somewhat larger overlap of 57%, but is based 
on a very different wave function and technique, i. e., a 
comparison of deformed-basis AMD with generator- 
coordinate method to a resonating-group wave function 
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0.25 




c. m. distance [fm] 



FIG. 7: Full curve: overlap of the strongly deformed state 
found in 32 S with a configuration of two 16 O nuclei placed 
symmetrically at a given distance. Dashed curve: norm of 
the combined Slater determinant of the two 16 O nuclei. 



for the double- 16 system. It is not surprising that a 
mean-field description of this strongly deformed state 
contains less of an na-cluster configuration than the 
AMD description. The AMD wave function, on the other 
hand, was found to have a larger overlap with the mean- 
field one, also its spin-orbit energy is in surprisingly good 
agreement with the mean-field one. This could be an in- 
teresting point for future studies. 

A similar attempt to increase the overlap with the na- 
clustcr configuration by varying the distance between the 
two quadruplet positions brought no further improve- 
ment. 



F. The Compound System 36 Ar 

This is somewhat similar in showing a slight necking-in 
effect, and again the iterations procedding from the two 
different initial configurations 16 O+ 20 Ne and 12 C+ 24 Mg 
led to the same final state. The shape is shown in Fig. [8] 
It is interesting that for this system the calculation with 
pairing found this superdeformed state rapidly from a 
wide range of initial distances, whereas the unpaired cal- 
culation tended to get stuck near j3 = 1 and only relaxed 
to the correct state after several thousand iterations. 

The unconstrained a-cluster analysis showed a diffuse 
picture. The configuration yielding the overlap quoted 
had two particles at ±2.7 fm, a quadruplet in a distorted 
rectangle at -0.3 fm and a triplet at +0.4 fm. The overlap 
depending only weakly on the precise positioning, this 
is also consistent, e. g., with the central seven particles 
placed in two triangles facing each other and a single a 
between at the center-of-mass. 




FIG. 8: Density contour lines for 36 Ar (left) and 40 Ca (right) 
superdeformed states. 



G. The Compound System 4U Ca 

This is really hyperdeformed as seen in Fig. |H1 The 
density shows a triple-humped structure, which is also 
borne out by the unconstrained placement of a-particles. 
It produced a configuration characterized by two parti- 
cles near ±4.3 fm, two particles at the center of mass, 
and two triangular triplets near ±2.9 fm along the sym- 
metry axis. The overlap was quite small, however, and 
seemed to depend only weakly on the exact placement 
of the particles, so that the individual particles were up 
to half a fermi from symmetric positions along the axis; 
the lateral displacement was always < 0.2 fm, though. A 
further indication that at this point the cluster interpre- 
tation becomes quite fuzzy is that for Sly6, which showed 
a practically identical density distribution, one of the a's 
from the triangles move closer to those on the end posi- 
tions. 



VII. EXOTIC a- PARTICLE STRUCTURES 

It should be briefly mentioned that the present code 
also allows the search for more exotic a-clustering struc- 
tures, such as chains and polygons, by initializing the 
static Mean — Field calculation with such a geometric ar- 
rangement. No interesting state was found for 3 and 4 
cv-particles, however, and as an example we confirmed the 
results of Itagaki and collaborators [26| that a three-a lin- 
ear chain is unstable with respect to triaxial deformations 
also for the mean field. 



VIII. SUMMARY 

The results obtained in this paper give a large amount 
of information on the relationship between cluster wave 
functions and ground-state Mean-Field ones. Let us sum- 
marize the main features of the ground-state analysis 
first: 

1. The automatic search for a na-cluster configura- 
tions with maximal overlap to the Mean-Field state 
worked surprisingly well and appears to really pro- 
duce the best configuration, as all analyses with 
specific geometric assumptions showed. The ge- 
ometry of the clusters appears similar to what is 
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produced in cluster calculations. 

2. Even when it was freely varied, the radius param- 
eter for the a's always has reasonable values. It 
slowly increases with the mass of the nucleus, indi- 
cating the diminishing localization of the particles. 

3. In all cases, the overlap decreased as the spin-orbit 
energy produced by a certain force parametrization 
increased. As expected the spin-orbit coupling de- 
stroys the spin-symmetric cluster structures. 

4. The overlap with clusters also decreased when go- 
ing to heavier systems. The growing effects of spin- 
orbit coupling and Coulomb explain that naturally. 

5. Also as expected the cluster wave functions worked 
better for the lighter systems up to 20 Ne. There 
was, however, again enhanced overlap for 28 Si and 
36 Ar. It is interesting that the AMD model deviates 
from a pure na-structure in the same way and stays 
closer to the mean-field model. 

6. In many cases the optimal na-structure contains 
groups of a-particles spaced very close together. 
They do not represent a true a-cluster configura- 
tion, but correspond to heavier clusters. 

7. The attempts at including correlations in the 
form of collective superpositions were disappoint- 



ing, producing only minor improvements in the 
overlaps. Apparently their relatively large size 
compared to the total nucleus enables the clusters 
to fill the nucleus without additional spreading. 

8. The relativistic model results were generally very 
similar to the Skyrme-force ones. The only remark- 
able difference is in the cluster radii, which were 
systematically smaller for the relativistic version. 

For the superdeformed states the clearest example of 
cluster structure is 32 S, and it is especially interesting 
that the overlap with the double- 16 system agrees very 
well with that of the unconstrained a-configuration. 
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